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ABSTRACT 


This  report  describes  an  investigation  of  the  effect  of  low  frequency 
periodic  perturbations  of  the  height  of  a  conducting  surface  on  the  earth's 
local  electric  field. 

Simulation  of  small,  low  audio  frequency  vibrations  of  a  portion  of  a 
conducting  surface  indicates  that  such  displacement  will  result  in  pertur¬ 
bations  to  the  normal  component  of  the  electric  field  above  that  surface. 

The  exact  magnitude  of  these  perturbations  is  a  function  of  the  type  of 
probe  used  to  sense  the  perturbations  and  the  probe-to-surface  distance 
and  is,  in  general,  on  the  order  of  microvolts. 

Atmospheric  noise  spectra  were  obtained  by  use  of  a  passive  antenna 
measurement  scheme.  The  data  obtained  were  analyzed  using  the  Fast  Fourier 
Transform  and  noise  spectra  on  the  band  0.1  to  100  Hertz  were  acquired.  The 
noise  spectra  thus  obtained  indicate  that  atmospheric  noise  appears  to  fall 
off  at  the  rate  of  27  dB  per  decade  in  the  region  0.1  to  100  Hertz. 

Examination  of  the  signal- to-noise  ratios  implies  that  the  performance 
of  this  system  in  detection  of  sea  surface  vibrations  would  be  marginal  at 
best . 
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INTRODUCTION 


In  the  past,  work  or,  Contract  Nonr-991  108 t  has  beer,  directed  toward 
various  techniques  for  sensing  „. .dewater  acoustic  signals  from  a  platform 
entirely  removed  from  the  surface  of  the  sea.  The  early  phases  of  the  work 
were  concerned  with  detecting  the  perturbations,  caused  by  acoustic  energy, 
to  the  permittivity  of  the  sea  water.  In  the  later  _.ears  :f  the  research, 
work  was  directed  towards  a  technique  to  detect  small  scale  vibrations  of  tie 
sea  surface,  again  caused  by  acoustic  energy  propagating  beneath  the  surface. 

In  this  phase  of  Contract  Nonr-991  f 08 )  ,  research,  has  been  conducted 
to  investigate  che  detectability  of  manmade  free  water  surface  motion  bv 
sensing  means  of  perturbations  to  the  normal  component  of  the  earth's 
atmospheric  electric  field. 

The  research  work  has  been  conducted  with  two  objectives  in  mind:  first, 
to  pcovide  a  model  for  the  sea  surface  motion  in  order  to  estimate  the  magni¬ 
tude  of  the  signal  one  could  expect  to  detecr,  and  second,  to  characterize 
spectrally  the  random  variations  in  the  earth's  atmospheric  electric  field 
tj^i  i  /*}■»  v'ould  3  s  b^ckr*!TO,_!!id  noiss  ir*.  ths  dst.se  tic  i'.  c  i  p  rr.  s  r.  t  "lb  see 

data  are  then  used  to  estimate  a  signal - to-noise  ratio  chat  characterizes 
the  detection  process. 

Since  the  tvpes  of  acoustic  activity  for  which  detection  is  desired  arc- 
expected  to  induce  sea  surface  vibrations  with  frequencies  general Iv  less 
than  100  Hertz,  the  noise  spectra  were  measured  in  the  range  of  O.i  to  100  Hertz. 


Section  II  of  this  report  discusses  the  estimation  of  vibration-induced 
signals  and  a  typical  detection  scheme.  Section  111  describes  the  noise 
characterization  technique,  while  Section  IV  analyzes  the  results  of  measure¬ 
ment  of  atmospheric  electric  field  fluctuations.  The  Appendices  describe 
specific  features  of  the  data  collection  and  analysis  procedures. 


II. 


SIGNAL  ESTIMATES 


A.  Background 

When  dealing  with  electrostatic  boundary-value  problems,  certain 
conditions  must  be  satisfied.  Some  of  these  conditions  are: 

(l)  Laplace's  equation  for  the  potential,  V,  must  equal  zero  every¬ 
where  except  on  boundaries, 

ine  tangential  component  of  the  elficn'ic  lield,  E  ^ ,  must  be 
continuous  across  the  boundary. 

(3)  The  normal  component  of  the  electric  flux  density,  D^,  must 

be  discontinuous  by  the  amount  of  surface  charge  density  at  the 
boundary. 

If  we  have  a  conducting  surface  with  a  periodic  perturbation  on  ft 
described  by 

Z  »  C  cos  £X 
o  "o 

and  static  electric  field 


condition  (2)  requires  that  the  field  be  perpendicular  to  the  surface. 
Solving  Laplace's  equation  for  the  potential 
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subject  to  the  boundary  conditions: 


V  -  0  when  Z  *  Z  -  f  cot  Bx 
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V  «  E  7.  when  Z  —  , 

o  ’ 


we  find  the  solution 
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When  P,'  «  1,  this  expression  can  be  simplified  to 
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Both  the  electric  field  and  the  potential  are  disturbed  in  the  same 
manner  as  Che  surface.  The  effect  of  the  surface  disturbance  falls  off 
exponentially  as  we  move  up  from  the  surface.  Therefore,  either  field 
measurements  or  potential  measurements  made  near  the  surface  will  show  the 
surface  perturbation.  Of  these  two  measurements,  the  potential  measurement 
is  the  easiest  tc  make  because  it  is  a  point  measurement,  while  the  E  field 
is  a  difference  or  gradient  measurement. 

One  way  that  these  measurements  can  be  made  is  to  suspend  a  sphere  oi 
a  conductor  of  any  other  practical  shape  above  the  Surface.  If,  in  the 
absence  of  the  sphere,  the  potential  of  this  point  is  V,  then  when  the 
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sphere  is  grounded  in  this  position,  the  conductor  must  transfer  a  charge 
of  such  value  as  to  make  the  potential  zero.  If  proper  instrumentation  Is 
placed  ir.  the  ground- to-probo  path,  the  charge  can  be  measured  and  the 
potential  determined.  If  the  change  in  potential  due  to  rnc  surface 
perturbation  is  all  that  is  desired,  then  the  distortion  of  the  field  caused 
by  the  probe  u i 1 1  not  have  to  be  accounted  for. 

This  type  of  measurement  can  also  be  used  to  measure  surface  displace¬ 
ments  whose  periods  are  shorter  than  the  relaxation  time  tor  electrical 
phenomena  in  the  atmosphere.  This  relaxation  time  is  defined  as  the  time 
it  takes  to  transfer  charge  tc  or  from  a  conductor  by  ionic  conduction. 

Since  actual  values  of  the  relaxation  time  for  air  near  the  earth's  surface 
are  from  5  to  40  minutes  depending  on  the  amount  of  pollution,^  the  case 
of  primary  interest  in  this  scudy  is  the  one  in  which  the  period  of  the 
displacement  is  shorter  than  the  relaxation  time.  In  this  case,  the  problem 
is  dealt  with  in  terms  of  electrostatic  effects  rather  than  current  effects 
and  the  change  of  potential  between  the  probe  and  ground  is  related  to  the 
change  of  the  probe-ground  capacitance. 

A  model  of  the  probe  suspended  above  a  conducting  plane  is  shown  m 
Figure  1.  In  this  figure  the  capacitor  is  the  probe-ground  capacitance, 
the  resistor  is  the  resistance  of  the  air  column  between  the  ground  and  the 
probe,  and  the  current  source  represents  the  conduction  current  that,  flows 
in  the  air.  This  model  is  valid  for  the  two  conditions  mentioned  above, 

If  a  perturbation  causes  the  probe-to-ground  distance  to  change  slowly  us 
compared  to  the  RC  time  constant,  then  the  capacitive  reactance  is  large 
compared  to  the  resistance  and  the  voltage  from  the  probe-to-ground  is  just 
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Figure  1.  Siepended  Probe  Model 
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the  potential  of  that  point  in  space.  However,  if  the  perturbation  causes 
the  probe- to-ground  distance  to  change  rapidly  with  respect  to  the  time 
constant,  the  capacitive  reactance  is  small  compared  to  the  resistance, 
the  probe-voltage  remains  constant;  and  we  have  a  electrostatic  problem. 


B.  Simulation  Model 

In  order  to  estimate  the  magnitude  of  the  field  disturbances  caused  by 
surface  height  perturbations,  a  simulation  model  of  the  proposed  system  was 
devised.  The  model  shown  in  Figure  1,  enlarged  to  include  a  few  more 
variables,  was  used  for  this  purpose.  A  schematic  diagram  of  this  circuit 
is  shown  in  Figure  2.  This  circuit  uses  a  Thevenin  equivalent  voltage 
source  in  place  of  the  current  source  shown  in  Figure  1,  and  the  capacitance, 
Cj,  and  resistance,  ,  represent  the  input  impedance  of  the  measurement 
instrumentation  and  the  stray  capacitance  of  the  leads.  Since  we  are 
primarily  interested  in  events  that  occur  faster  than  the  relaxation  time, 
the  potential  at  the  probe  is  shown  as  a  battery. 

Let  be  the  voltage  across  and  be  the  voltage  across  C^,  C ^ 
and  R£.  Then  solving  for  V wc  know  that 


X1  +  i2  +  *3 


dql  dq2  V2 

where  ^  >  i2  =  “  ’  13  “ 


ql  =  C1V2  ’  and  q2  ~  C2V? 
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If  the  surface  below  the  probe  Is  perturbed  in  some  periodic  manner. 
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where  is  the  capacitance  at  maximum  spacing 

b  is  the  ratio  of  peak-to-peak  displacement  to 
maximum  plate  separation,  and 

a'  is  the  angular  frequency  of  the  perturbation 

Equation  (5)  can  now  be  rewritten 


C.  Simulation  Results 

An  estimate  for  the  atmospheric  potential,  V,  can  be  derived  from  the 

fact  that  the  average  fair-weather  value  of  the  earth's  electric  field  is 

100  volts  per  meter.'  The  capacitance,  C^,  can  be  calculated  from  probe 

size  and  height  figures,  and  can  be  estimated  from  the  fact  that  the 

time  constant,  C^,  will  be  in  the  range  of  5  to  40  minutes,  depending 

upon  the  amount  of  pollution  in  the  air. 

For  a  1  meter  diameter  plate  mounted  15  centimeters  above  the  ground, 

the  probe- to-grcund  capacitance  would  be  approximately  50  picofarads  and  the 

potential  between  the  plate  and  the  ground  would  be  approximately  15  volts. 

14 

Using  these  values  for  Cq  and  ,  and  C 2  =  122  picofarads  and  *  10  ohms, 
the  input  impedance  and  cable  capacitance  values,  we  can  solve  equation  (7). 

Figure  3  shows  the  response  as  a  function  of  displacement  for  fixed 
values  of  initial  capacitance  and  Figure  4  shows  the  voltage  response  as 
a  function  of  initial  capacitance  for  given  values  of  displacement.  Exami¬ 
nation  of  these  figures  shows  that  the  output  falls  off  almost  linearly 
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Figure  3.  Probe  Response  as  Function  of  Di  ’aceaent,  f,  for 
Given  Values  of  Initial  Capacitance  'picofarads) 
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with  displacement  and  is  a  linear  function  of  capacitance  for  mid-values, 
while  it  is  a  higher  order  function  at  the  extremes.  We  can  see  from 
Figure  3  that,  for  displacements  on  the  order  of  a  few  thousandths  of  a 
percent  of  the  probe-to-surface  spacing,  the  voltage  change  is  going  tc 
be  in  the  microvolt  region. 

This  model  can  be  further  refined  by  having  only  a  portion  of  the 
bottom  plate  vibrate.  This  resembles  the  real  world  case  more  than  the 
previous  model  where  the  whole  bottom  plate  vibrated. 

Figure  5(a)  shows  the  geonetry  cf  the  capacitor  with  pert  of  the 
bottom  plate  at  maximum  displacement.  This  configuration  can  be  trans¬ 
formed  to  the  normal  parallel  plate  capacitor  shown  in  Figure  5(b)  through 
the  use  of  a  Schwarz-Christof fel  transformation  (see  Appendix  A).  Using 
this  transformation  assumes  that  the  plates  are  infinite  planes.  However, 
if  the  spacing  between  the  plates  is  small  compared  to  the  width  of  the 
plates,  then  the  solutions  obtained  are  good  approximations  of  the  finite 
case.  We  can  extend  our  approximation  to  the  case  shown  in  Figure  5(c) 
by  taking  a  truncated  mirror  image  cf  Figure  5(a)  and  adding  it  to  the 
original  capacitor.  This  will  be  valid  so  long  as  h  is  small  compared  to  d. 

Transforming  a  stepped-capacitor  to  the  w-plane  and  solving  for  the 
capacitance,  we  find  that  the  increase  in  capacitance  is  slightly  more  than 
that  obtained  from  two  parallel  capacitors,  one  at  the  wide  spacing  and  one 
at  the  narrow  spacing.  This  is  expected  because  the  transformation  takes 
the  internal  fringing  into  account  which  is  neglected  when  we  consider  the 
stepped  capacitor  as  two  parallel  ones.  If  we  further  consider  the  stepped- 
capacitor  shown  in  Figure  5(c)  as  the  cross  section  of  three  disks  instead 
of  infinite  planes,  the  error  in  capacitance  value  is  greater  due  to  the 

i 


13 


increased  amount  of  fringing  that  is  being  neglected.  However,  considering 
the  capacitor  as  a  fixed  one  and  a  smaller  variable  one  will  give  us  an  order 
of  magnitude  results. 

For  example,  assume  we  have  a  one-meter  diameter  probe  at  a  height  of 

15  centimeters  above  the  surface;  the  capacitance  of  this  probe  will  be 

approximately  50  picofarads.  If  the  area  of  perturbation  has  a  radius  of 

5  centimeters  and  a  peak-to-peak  displacement  of  1.5  millimeters,  this  will 

give  us  an  initial  capacitance  of  approximately  0.5  picofarads  with  5  equal 

to  0.01.  Referring  to  either  Figure  3  or  4  with  these  values  we  find  that 

-4 

the  voltage  perturbation  will  be  on  the  order  of  10  volts. 

A  laboratory  model  of  the  parallel  plate  vibrating  capacitor  was 
constructed  to  test  the  validity  of  the  measurement  procedure.  The  model 
consisted  of  two  parallel  plates  6  inches  in  diameter  connected  to  a  voltage 
source  through  a  resistor.  The  bottom  plate  was  attached  to  a  cam  and  driven 
by  a  motor  at  a  30  Hz  vibration  frequency.  The  experimental  results  for 
large  values  of  6  agreed  with  the  theoretical  value  of  within  10-207.  and 
to  within  1 07  for  small  values  of  5.  Most  of  this  error  was  attributed  to 
mechanical  instability  of  the  modsl  and  qualitative  experiments  were  run  to 
prove  that  the  varying  component  of  the  voltage  was  in  fact  caused  by  the 
varying  capacitance  and  not  by  any  outside  influences. 
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HI.  NOISE  ESTIMATES 


A.  Background 

In  Section  II-A  we  showed  how  a  surface  disturbance  could  cause  the 
potential  at  some  point  above  it  to  vary.  Me  further  described  a  system 
to  measure  these  disturbances.  As  in  all  physical  systems,  there  are 
exogenous  factors  that  enter  the  picture  as  noise  or  unwanted  influences. 
The  factors  to  be  considered  in  this  system  are  man-made  electrical  noise, 
such  as  60  Hz  noise,  and  random  fluctuations  in  the  local  earth's  electric 
field. 

Using  Maxwell's  equations  and  assuming  that  the  electric  field  Is 
static  and  vertical  and  depends  only  on  altitude,  we  know  that 


J  ■  aE 


and 


dj 

dz 


therefore 


j  -  j  “  constant 
o 


where  o  is  the  charge  density,  e  the  free  space  permittivity,  a  the  conduc- 
tivity,  E  the  electric  field  strength,  and  j  the  current  density.  It  can  be 
shown  that,  subject  to  the  boundary  conditions 

z  -  od  ,  V  -  0, 
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the  potential  at  the  surface  of  the  earth  is 

a  q 

V(0)  *  —  R 

e 

o 

where  q  is  the  charge  on  the  earth  and  R  is  the  total  resistance  of  the 
atmosphere.  This  equation  has  the  form 


V  -  IR 

where 


Variations  in  this  current  can  be  caused  by  changes  in  the  conductivity  or 
the  surface  charge  density.  The  conductivity  varies  with  the  amount  of 
pollution  in  the  air  and  with  changes  in  ionization  levels  due  to  increased 
cosmic  ray  penetration,  point  discharges,  a.'.d  lightning  flashes.  The  surface 
charge  density  can  vary  locally  with  changing  upper  atmosphere  charge  con¬ 
ditions  such  as  cloud  passage  and  precipitation. 

In  order  to  examine  the  noise  mechanism,  the  simplified  model  used 
the  signal  analysis  has  to  be  expanded.  From  the  above  model  we  can  see  that 
the  current  source  is  connected  to  the  probe  through  a  parallel  resistor- 
capacitor  network.  The  expanded  model  is  shown  in  Figure  6.  This  model  fits 
the  actual  atmospheric  situation  quite  well.  Suppose  the  lower  atmosphere  is 
in  a  quasi-static  state  and  a  fresh  charge  cloud  suddenly  appears.  This  will 
immediately  produce  a  potential  gradient  at  the  ground  (or  any  t  ther  reference 
point),  calculable  purely  electrostatically .  But,  also  immediately,  the 
current  flow  will  be  altered,  charges  will  move  and  a  new  quasi-static  state 
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will  be  set  up.  This  situation  can  be  illustrated  using  Figure  6. 
change  in  charge,  Aq,  causes  a  corresponding  change  in  voltage,  AV^  and 
AV-,  across  the  capacitors.  Since  these  are  parallel  RC  circuits,  the 
voltages  across  them  at  any  time,  t,  due  to  this  disturbance  are  equal  to 

-t/R  C 

v|(t)  =  AVj  e 

-t/R  CL 

Vj ( t)  -  av2  e 


The  change  in  charge  alsc  changes  the  value  of  the  current.  This  produces 
a  voltage  equal  to 


V”(t) 

V’^(t) 


AIR^l  - 
AIR2(1  - 


-t/R.C. 

4  i) 


-t/R_C, 
*  Z  ) 


Therefore  the  voltage  at  the  probe  at  any  given  time  is  the  sum  of  these  or 


v  (t)  =  v;  (t)  +  vM(t) 

p  i  i 

-t/R.C  -t/R  C 

V  (t)  *  AV,  «  1  1  +  AIR, (1  -  * 

pi  1 

If  there  is  a  change  of  conductivity  rather  than  of  charge,  then  there  is  no 
electrostatic  effect  and 
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Further  examination  of  this  model  shows  that  the  probe  potential  follows 


that  of  the  surrounding  atmosphere  quite  rapidly.  If  we  take  the  special 
case  where  the  conductivity  is  constant  (i.e.,  no  space  charge  distribution), 
this  circuit  acts  as  a  simple  resl3tive  divider  network  since  ■  ^2^2’ 

Furthermore,  if  the  probe  is  a  thin  disk  parallel  to  the  earth's  surface,  then 
the  field  distortion  due  to  it  will  be  confined  to  the  immediate  proximity  of 
the  probe.  We  can  then  consider  the  upper  charge  layer  to  lie  on  any  undis¬ 
turbed  equipotential  line  close  to  the  probe.  This  being  the  case,  even  with 
conductivity  varying  with  height,  it  is  readily  apparent  that,  with  all  the 
changes  assumed  to  be  voltage  changes  occurring  at  some  undisturbed  equi¬ 
potential  line  close  to  the  probe,  the  probe  will  Immediately  respond. 


B.  Experiment  PesiRn 

To  measure  and  record  the  background  fluctuations  in  the  vertical  com¬ 
ponent  of  the  earth' 6  electric  field,  a  probe  consisting  of  a  one-meter 
diameter  plate  was  mounted  15  centimeters  above  a  ground  plane  approximately 
four  meters  square.  The  plate  diameter  and  spacing  were  selected  to  provide 

a  capacitance  of  approximately  50  picofarads  and  an  equivalent  leakage 

12 

impedance  of  approximately  10  ohms.  With  this  leakage  impedance,  the 

14 

input  impedance  of  the  electrometer  (po  10  ohms)  is  large  compared  to  the 
leakage  impedance  and  thus  the  effect  of  the  electrometer  may  be  neglected. 


Figure  i  shows  a  pictorial  diagram  of  this.  In  this  drawing,  Cq  represents 
the  capacitance  between  the  probe  and  ground,  represents  the  input  capacitance 


ot  the  electrometer  and  the  capacitance  of  the  cables  connecting  the  probe 
to  the  electrometet ,  and  R^  represents  the  input  resistance  of  the  electro¬ 
meter. 

Several  schemes  were  investigated  for  coupling  the  electrometer  to  the 
probe.  The  method  which  proved  most  successful  used  a  triaxial  cable  made 
from  a  length  of  low  loss  coaxial  cable  with  an  additional  outer  shield. 

The  inner  shield  was  used  as  a  guard  in  conjunction  with  a  feedback  circuit 
to  reduce  the  effect  of  the  cable  capacitance.  The  circuit  used  with  the 
electrometer  is  shown  in  Figure  8,  As  can  be  seen,  an  output  signal  from  the 
electrometer  proportional  to  the  voltage  between  the  center  conductor  and 
guard  is  amplified  in  the  operational  amplifier  and  fed  back  to  the  guard, 
thus  reducing  the  difference  in  potential  between  the  center  conductor  to 
a  level  determined  by  the  gain  of  electrometer  and  the  operational  amplifier 
circuitry.  This  difference  in  potential,  with  the  equipment  used  in  this 
experiment,  was  a  maximum  of  30  millivolts  with  the  center  conductor  (which 
is  connected  to  the  probe)  at  120  volts  above  ground.  Since  the  difference 
in  potential  between  the  center  conductor  and  guard  is  small  compared  to  the 
difference  in  potential  between  the  central  conductor  and  ground,  we  can 
measure  the  potential  difference  between  the  guard  and  ground  and  use  this 
as  c.n  approximation  to  the  actual  potential  difference  on  the  probe  with  an 
error  of  less  than  0.037*.  In  addition,  the  low  output  impedance  (current 
source)  characteristics  of  the  operational  amplifier  make  it  useful  for 
driving  the  recording  instruments.  An  analysis  of  the  gain  of  the 
electrometer-feedback  circuit  is  included  in  Appendix  B.  An  important 
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feature  of  the  feedback  circuit  is  the  reduction  of  the  effective 
capacitance  of  the  cable.  As  is  shown  in  the  analysis  in  Appendix  C,  the 
effective  capacitance  is  reduced  by  approximately  the  gain  of  the  electrometer- 
feedback  circuicrv.  The  effective  input  capactiance  of  the  electrometer  and 
cable  to  small  signal  changes  is  approximately  0.5  picofarad  which  is  two 
orders  of  magnitude  less  than  the  capacitance  of  the  probe,  and  therefore 
has  an  inconsequential  effect  on  the  measurements, 

'he  frequency  response  of  the  probe  can  be  determined  by  referring 
back  to  the  model  in  Figure  6.  For  time  variations  slow  compared  to  the 
time  constant  of  the  atmosphere,  the  transfer  function  for  this  circuit  will 
be  essentially  that  of  a  resistive  voltage  divider  and  is  frequency  independcr, t . 
For  time  variations  on  the  order  of  the  relaxation  time,  the  transfer  func¬ 
tion  will  be  that  of  the  series  parallel  RC  circuit  and  a  strong  function 
of  frequency.  In  the  third  regime,  however,  time  variations  are  much  faster 
than  the  atmospheric  time  constant;  thus  the  capacitive  reactance  will  be 
small  compared  to  the  resistive  component.  The  circuit  will  behave  as  a 
capacitive  voltage  divider,  and  the  transfer  function  will  again  be  essentially 
independent  of  frequency.  Thus,  if  frequencies  much  greater  than  the  reciprocal 
of  the  atmospheric  relaxation  time  are  of  interest,  no  consideration  need  be 
made  of  the  frequency  response  of  the  probe.  The  primary  considerations  in 
probe  design  are  that  the  input  impedance  of  the  electrometer  be  high  com¬ 
pared  to  the  leakage  resistance  at  the  probe  and  that  the  ratio  of  probe 
capacitance  to  miscellaneous  stray  capacitance  (such  as  cable  capacitance) 
be  made  large  as  feasible.  This  is  due  to  the  fact  that  cable  capacitance 
increases  the  apparent  sise  of  the  probe-to-ground  capacitance  without  allowing 


additional  coupling  to  the  atmospheric  electric  field  which  reduces  the 
usable  sensitivity  of  the  probe  to  variations  in  the  atmospheric  electric 
field. 

In  order  to  obtain  the  dynamic  range  needed  to  display  the  whole 
spectrum  from  0  to  100  Hz,  some  signal  pre-conditioning  was  needed.  This 
was  accomplished  through  the  use  of  filters  and  amplifiers  to  selectively 
amplify  portions  of  the  spectrum.  Active  filters  of  two  types  were  con¬ 
structed.  One,  a  six  pole  Butterworth  low  pass  with  a  100  Hertz  break 
frequency  was  used  to  band  limit  Che  data  for  purposes  of  eliminating 
aliasing  in  the  analysis  program.  A  second  filter  was  constructed,  a  0.01 
Hertz  six  pole  active  Chebyshev  highpass,  and  followed  by  an  amplifier 
stage  having  a  10  dB  gain.  In  recording  the  data,  two  channels  on  the 
recorder  were  used;  one  channel  recorded  the  low-passed  data,  while  the  low- 
passed  data  passed  through  the  0,01  Hertz  high  pass  filter  and  amplifier 
was  recorded  on  another  channel.  Data  records  of  thirty  minutes  to  two 
hours'  duration  were  made.  After  recording,  the  data  was  converted  to  digital 
form  on  a  Radiation  Analog-to-Dlgital  Converter.  The  equivalent  sampling 
rate  used  in  A/D  conversion  was  1000  samples/second,  thus  insuring  that 
all  signal  components  to  100  Hz  were  recovered  unambiguously.  After  the 
A/D  operation,  the  digitized  data  was  spectrally  analyzed  by  Fast  Fourier 
Transform  techniques. 

Figure  9  shows  some  typical  results  obtained  from  th< se  measurements. 

These  were  obtained  on  a  clear  day  with  a  slight  wind.  Figure  10  shows  some 
results  obtained  when  there  was  thunderstorm  activity  in  the  vicinity,  the 
sky  was  cloudy  and  there  was  a  slight  wind.  These  results  are  somewhat  typical 
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Figure  10.  Atmospheric  Potential  Fluctuations  as  Measured  by  1  Meter 

Diameter  Probe  15  cm  above  Ground  Plane — Thunderstorm  Activity 
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and  are  shown 


12  3 

and  appear  to  agree  with  other  published  records,  ’  ’ 
principally  to  demonstrate  that  the  system  is  capable  of  following  rapid 
changes  and  is  stable  during  slow  variations. 

C.  Experimental  Results 

There  are  several  options  available  for  analysis  of  the  data.  One 
is  to  calculate  the  spectrum  via  the  correlation  function  of  the  data;  a 
second  way  is  to  directly  transform  the  data.  With  either  of  these  two 
approaches,  weighting  functions  (lag  windows  and  data  windows)  can  be  used 
to  modify  the  data. 

The  first  method  above  is  the  classical  method  of  computing  power 
spectra.  One  of  the  reasons  for  its  widespread  use  was  its  computational 
economy.  It  took  markedly  fewer  operations  to  calculate  perhaps  a  tenth  as 
many  mean  lagged  products  as  there  were  data  points  and  to  Fourier  transform 
the  results  than  it  did  to  calculate  all  the  Fourier  coefficients.  The  advent 
of  the  Fast  Fourier  Transform  (FFT)  has  changed  that — it  is  now  faster  to 
calculate  all  the  Fourier  coefficients.  The  former  method  is  still  highly 
popular,  however,  because  it  is  so  well  established  and  most  detailed 
practical  problems  encountered  have  accepted  solutions  available. 

When  all  the  Fourier  coefficients  are  calculated,  it  would  appear  that 
the  best  approach  would  be  to  calculate  the  raw  periodogram.  This  ia  basically 
a  good  approach,  particularly  if  we  have  a  periodic  signal.  However,  with 
noise  type  signals,  the  leakage  involved  with  the  finite  Fourier  transform 
is  too  high.  The  best  way  to  attack  this  problem  is  through  the  use  of 
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windows.  The  choice  of  a  window  depends  on  the  bandwidth  of  the  effective 
filter  and  the  statistical  stability  desired  ir  the  estimate.  In  other 
words,  we  are  concerned  with  the  trade-off  between  leakage  and  frequency 
resolution.  The  raw  periodogram  gives  the  best  resolution  and  the  worst 
leakage.  Window  functions  reduce  the  leakage  but  at  the  expense  of 
resolution.  This  is  due  to  the  fact  that  multiplying  the  data  by  a  data 
window  is  equivalent  to  spectral  averaging  in  the  frequency  domain.  The 
mean  square  error  is  reduced  but  the  spectral  estimate  is  smeared  out. 

Several  window  functions  were  considered  and  tested.  One  of  the 
easiest  to  use,  and  also  one  of  the  most  popular,  is  cosine  bell  data 
window  or  Hann  window.  This  corresponds  to  a  three  point  smoothing  with 
the  points  weighted  1/4,  1/2,  1/4.  The  effective  bandwidth  of  the  Hann 
window  is  twice  that  of  the  raw  periodograra  and  the  leakage  is  considerably 
reduced  due  to  the  reduction  of  the  side  lobes. 

Another  popular  window  function  is  the  Hamming  window.  This  corres¬ 
ponds  to  three  point  smoothing  with  weights  0.23,  0.54,  0.23.  The  bandwidth 
is  just  about  the  same  as  the  Hann  window,  the  first  side  lobe  is  lower  but 
the  height  of  the  aide  lobes  do  not  drop  as  fast  as  they  do  with  the  Hann 
window. 

A  third  window  function  is  the  Parzen  window.  The  bandwidth  of  this 
function  is  about  25%  greater  than  that  of  the  Hann  or  Hamming  windows 
and  the  side  lobes  are  substantially  reduced.  Because  the  main  lobe  is 
wider,  the  variability  of  the  corresponding  spectral  estimates  is  less. 

Since  the  use  of  window  functions  is  equivalent  to  smoothing  the  frequency 


29 


domain  data,  five  point  smoothing  with  equal  weights  was  also  used  in  the 
analysis.  The  equivalent  bandwidth  in  this  case  is  five  times  that  of  the 
raw  periodogram;  this  reduces  the  variability  of  the  spectral  estimates 
considerably.  One  measure  of  this  variability  is  the  mean  square  error  and 
it  was  calculated  for  each  of  the  windows  and  found  to  be: 


e 

(raw  periodogram)  = 

1.0 

e 

(Hann  window)  • 

0.707 

c 

(Hamming  window)  = 

0.707 

c 

(5  pt.  smoothing)  = 

0.447 

e 

(Parzen)  ■ 

0.577 

The  information  of  most  interest  to  this  study  was  the  rate  of  roll  off 
of  the  spectrum.  Therefore,  statistical  variability  of  the  estimates  was  of 
greater  importance  than  the  resolution  bandwidth.  For  this  reason,  the 
Parzen  window  and  the  five  point  smoothing  appeared  to  be  the  best  candidates 
for  window  functions.  The  data  was  analyzed  using  these  window  functions  and 
the  spectra  obtained  were  virtually  identical.  An  atmospheric  noise  spec¬ 
trum  obtained  from  these  measurements  is  shown  in  Figure  11.  The  system 
noise  was  also  measured  and  anlyzed  and  its  spectrum  is  shown  in  Figure  12. 

The  svstero  noise  data  were  obtained  with  the  probe  connected  but 

electrostatically  shielded  from  the  earth’s  electric  field. 

The  atmospheric  noise  spectrum  falls  off  at  a  rate  of  27  dB  per  decade 
from  0.2  Hz  to  100  Hz.  The  spectrum  level  remains  above  that  of  the  system  noise 
from  0.1  Hz  to  about  50  Hz  where  it  becomes  noise  limited.  The  major  problem 
seems  to  be  the  dynamic  range  available  on  the  tape  recorder  and  could  be 
overcome  by  further  selective  filtering  and  amplifying  of  the  data.  These 
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Figure  11.  Atmospheric  Nolae  Spectrum  from  0.1  to  100  Hertz 
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age 


results  are  adequate  to  complete  tue  analysis  of  the  proposed  surface 

perturbation  measurement  system  and  match  favorably  with  those  obtained 

A 

at  higher  frequencies  by  E.  L.  Maxwell.  His  results  show  the  field 
intensity  to  be  70  dB  below  1  volt/meter  at  100  Hs  and  falling  off  at 
about  20  dB  per  decade  above  that  point. 
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IV.  ANALYSIS  AND  CONCLUSIONS 


The  signal  estimates  and  the  nu.sc  es.timatts  car.  :.ow  I»e  compared  „nd 
some  conclusions  reached  about  the  feasibility  of  using  a  suspended  probe 
to  measure  surface  perturbations.  A  suitable  configuration  for  this  analvsis 
will  be  the  same  as  that  used  for  the  noise  measurements.  This,  it  will  be 
remembered,  consisted  of  a  1  meter  diameter  plate  suspended  15  centimeters 
above  the  ground,  and  the  capacitance  between  the  plate  and  g  ound  was  50 
picofarads,  tt  was  shown  in  Section  I  that,  with  this  arrangement,  a  pertur¬ 
bation  with  about  a  5  centime  ter  radius  and  a  peak-to-peak  displacement  cf 
1.5  millimeters  would  cause  a  voltage  perturbation  cn  the  order  of  10  **  volts. 
Referring  to  the  noise  spectrum,  we  see  that  at  10  Ha  the  noise  field  inten¬ 
sity  is  about  40  dB  below  1  volt/meter;  this  yields  a  noise  voltage  across 
the  15  centimeter  gap  of  about  10  ^  volts. 

With  this  highly  optimistic  situation,  the  sigoal-to-ncise  ratio  is  less 
than  1.  When  the  degrading  effects  of  the  real  world  are  taken  into  account 
in  the  signal  estimates,  the  ratio  will  be  even  worse.  Also,  most  practical 
situations  would  require  that  the  probe  be  suspended  1  meter  or  more  above 
the  ground.  This  would  reduce  the  capacitance  in  the  above  example  by  a 
factor  of  7  and  would  reduce  the  voltage  perturbation  by  two  orders  of  magni¬ 
tude.  The  effect  is  compounded  because  the  noise  voltage  is  increased  by  a 

-6 

factor  of  about  7.  We  now  have  a  signal  level  of  about  10  volts  and  a  noise 

-2 

level  of  about  10  volts.  Some  relief  can  be  obtained  by  going  to  a  higher 
frequency,  but  this  does  not  give  enough  because  the  noise  is  reduced  by  only 
an  order  of  magnitude  per  decade;  most  signals  of  interest  vould  be  below 


1000  Ha. 
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Appendix  A 


The  Schvarz-C’nristoffel  transformat ion  allows  the  mapping  of  fields 
with  polygonal  boundaries  in  the  o- plane  into  homogeneous  fields  with 
parallel  boundaries  in  the  v-plane.  In  most  cases  this  cannot  be 
accomplished  in  a  single  step,  and  it  is  first  necessary  to  map  the 
t-plane  on  an  auxiliary  t-plane  such  that  the  boundaries  of  the  z-plane 
go  over  into  the  real  axis  of  the  t-plane.  The  t-plane  is  then  transformed 
to  the  v-plane . 

The  problem  is  to  map  the  polygonal  boundary  of  the  o- plane  on  the 

real  axis  of  the  t-plane,  and  then  to  map  the  parallel  strip  of  the 

v-plane  on  the  t-plane.  Then  integration  constants  have  to  be  selected 

in  such  a  fashion  as  to  bring  about  a  direct  correspondence  between  these 

i  cy 

two  mappings.  If  we  let  (t  -  r^)  »  ♦  where  the  are  fixed 

points  on  the  real  axis  of  the  t-plane,  the  real  axis  of  the  t-plane  will 
transform  into  a  polygon  on  the  t-plane  since 


dz 

at 


.  (  ',-dn 

~.(c  -  d 


(t  -  r0) 


-ay 


does  not  change  directions  except  at  U.e  singular  points  r  ,  and  then  it 
changes  through  an  angle  a.n.  Expressed  in  integral  form  this  mapping 


function  is 


A  =  A 


1  J 


at 


(t  -  a1)ai  (t  -  a0)32 


Similarly,  the  t-plane  can  be  mapped  on  the  v-plane  by 


w 


_  dt 

(t  -  r/1  (t 


+  b2 
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figure  A-l  shows  the  boundaries  in  the  z- ,  t- ,  and  w-planes;  the 
Roman  numerals  indicate  the  points  where  two  sides  of  the  polygon  meet. 
The  unclosed  portions  are  assumed  to  close  at  infinity  and  are  indicated 
by  dotted  paths.  The  points  II  and  III  are  arbitrarily  located  in  the 
t-plane  at  0  and  +1  respectively.  The  parallel  boundaries  in  the  w-planc 
are  taken  to  be  equipotential  surfaces  with  v  =  0  and  v  =  V.  Setting  up 
the  differential  equations 


02 

dt 


(t  -  0)  (t  -  l)"2  (t  -  a)2 


r  1  /tTT  . 

2  =  J  ~  v  TTg  at  +  *2 


\  log  ^  . 
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Integrating  over  paths  I  to  determine  A^  and 
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Therefore,  -  JH  =  jA^n  for  R  -*  ® 


\  =  -« A 
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Upon  substitution  of  the  integration  constants 


AO 


APPENDIX  B 


The  actual  feedback  circuit  Is  constructed  around  a  Phllbrick  K2-XA 
vacuum  tube  operational  amplifier.  The  circuit  may  be  analysed  as  shown 
below. 


(3) 


(4) 
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Substituting  back  into  (4) 


(5) 


Thus 


^2  R3  ”  ^1 

_R  .  .  -jL(_2 - L  +  1} 

E2  R1  R4 


(6) 


The  output  of  the  feedback  circuit  is  feedback  along  the  low  side  of 
the  electrometer  input  to  form  a  guard  or  equipotential  surface  at  the  same 
(or  nearly  the  same)  potential  as  the  high  side  of  electrometer  input,  thus 
reducing  the  equivalent  input  capacitance.  The  output  of  the  feedback  circuit 
is  also  used  to  drive  the  recording  instruments.  Recordings  are  made  of  the 
output  voltage  on  a  standard  pen  chart  recorder  and  through  a  low  pass  filter 
on  magnetic  tape.  The  magnetic  tape  is  converted  to  digital  form  and  the 
spectra  computed  using  Fast  Fourier  Transform  techniques. 
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APPENDIX  C 


Ground 


Capacitance  of  cable 


CQ  center  to  guard 

Voltage  from  center  conductor  to  ground 
Voltage  from  center  conductor  to  guard 
Voltage  from  guard  to  ground 


from  Circuit  Analysis 
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If  we  understand  quasi-static  to  mean  that  there  are  to  be  negligible 
currents  flowing  from  the  probe  to  the  measuring  instrument;  or,  in  other 


Thus,  viewed  from  the  input  to  the  cable  (the  probe)  by  using  a 
guarded  cable  in  which  the  guard  is  maintained  at  nearly  the  same  potential 
(with  respect  to  ground)  as  the  center  conductor,  we  can  reduce  the  effect 
of  the  capacitance  of  the  cable. 


45 


Unclassified 
Security  Classification 


DOCUMENT  CONTROL  DATA  -  R&D 


cltHlftWlwi  »l  My  o*  •>•«<•«  and  mdai ^  «Mwt«Hon  omN  6a  inland  safewt  ft*  ewnll  i«p»N  I*  cRM*r(t*4> 


I  OaiGiwailwc  ACTIVITY  (Cerpcitr,  mithot) 

Georgia  Institute  of  Technology 

Atlanta,  Georgia  30332 

2 •  RCAOM  ItCuRltv  CLAIliflCATlON 

Unclassified 

2b  croup 

1  REPORT  TITLE 

INVESTIGATION  OF  MANMADE  AND  NATURAL  PERTURBATIONS  OF  THE  EARTH’S  SCIENTIFIC  FIELD 

*  DESCRIPTIVE  NOtES  fTyp«  of  report  and  tacluaiv*  c&a  re* ) 

Technical  Report  No.  6 

S  AUTHORfSj  f&bbt  ftrei  nm.  inttimJ) 

Bonham,  C.  H.  Ill;  and  Langley,  J.  B. 

ft  REPO  RT  DATE 

30  October  1970 

?■  TOTAL.  NO  OF  PASt*  ]  7  *>  NO  OF  RIM 

45  I  4 

•  «  CONTRACT  ON  (RANT  NO. 

Nonr-991 (08) 

b  RROJlC T  NO 

c. 

d 

9*  OR)  OiN  A  10R*>  Rt»ORT  N 

Technical  Report  No.  & 

A-5I9 

9b  OlHIR RtPOMT 
*»!•  report) 

_ 

Of (A  njr  otfrer  nsjer*6mrm  #vai  m>  b*  Aacifrtod 

10  AVAILABILITY,  LIMITATION  NOTICES 

11  SUPPLEMENTARY  NOTES 

k 

r 

- ^ - - 

SPOHSOSING  HUJTAHY  HCTIYITY 

Office  of  Naval  Research 

Department  of  the  Navy 

Arlington,  Virginia  22217 

15  ABsnyfer 


This  report  describes  an  investigation  of  the  effect  of  low  frequency 
periodic  perturbations  of  the  height  of  a  conducting  surface  on  the  earth's 
local  electric  field. 

Simulation  of  small,  low  audio  frequency  vibrations  of  a  portion  of  a 
conducting  surface  indicate  that  such  displacement  will  result  in  perturbations 
to  the  normal  component  of  the  electric  field  above  that  surface.  The  exact 
magnitude  of  these  perturbations  is  a  function  of  the  type  of  probe  used  to 
sense  the  perturbations  and  the  probe-.to-surface  distance  and  is,  in  general, 
on  che  order  of  microvolts. 

Atmospheric  noise  spectra  were  obtained  bv  use  of  a  passive  antenna 
measurement  scheme.  The  data  obtained  were  analysed  using  the  Fast  Fourier 
Transform  and  noise  spectra  on  the  band  0.1  to  100  Hercr  were  acquired.  The 
noise  spectra  thus  obtained  indicate  that  atmospheric  noise  appears  to  fall 
off  at  the  rate  of  27  dB  per  decade  in  the  region  0.1  to  100  Hert*. 

Examination  of  the  signal-to-noise  ratios  implies  thac  the  performance 
of  this  system  in  detection  of  sea  surface  vibrations  would  be  marginal 
at  best. 
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